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Abstract 

A generalized quantum distribution function is introduced. Tlie corre- 
sponding ordering rule for non-commuting operators is given in terms of 
a single parameter. The origin of this parameter is in the extended canoni- 
cal transformations that guarantees the equivalence of different distribution 
functions obtained by assuming appropriate values for this parameter. 



1 Introduction 

The hamiltonian formulation of classical mechanics treats the generalized coor- 
dinates, q, and momenta, p, as independent variables and in a symmetric way. 
In quantization procedure, however, one chooses the q- or p-representation de- 
stroying the equal standing of the two variables. In his pioneering work, Wigner 
[1] proposed quantum distribution function in phase space keeping the symme- 
try of q and p and defined the expectation values of observables in the manner 
of classical statistical mechanics. Here, we use the extended phase space tech- 
nique proposed by Sobouti and Nasiri [2] and the canonical transformations in 
this space to introduced a generalized quantum distribution function. It is known 
that for any distribution function there is an ordering rule of non-commuting op- 
erators. Thus, for generalized distribution function a generalized ordering rule 
is introduced, as well. Assuming appropriate values for the parameter involved, 
different distribution functions and their corresponding ordering rules can be 
obtained. The origin of these parameters is in the extended canonical transfor- 
mations in the extended phase space. 

The layout of the paper is as follows: In section 2, a brief review of the 
extended phase space formalism is presented. In section 3, the generalized dis- 
tribution function is introduced. Section 4, is devoted to conclusions. 
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2 A Review of EPS formalism 

A direct approach to quantum statistical mechanics is proposed by Sobouti and 
Nasiri [2], by extending the conventional phase space and by applying the canon- 
ical quantization procedure to the extended quantities in this space. Here, a brief 
review of this formalism is presented. For more details, the interested reader 
may consult Sobouti and Nasiri [2]. 

Let £'(5, q) be a lagrangian specifying a system in q space. A trajectory of 
the system in this space is obtained by solving the Euler-Lagrange equations for 
lit), 

1^_^=0. (1) 
dt dq dq 

The derivative calculated on an actual trajectory, that is, on a solution of 
Eq. (1), is the momentump conjugate to q . The same derivative calculated on a 
virtual orbit, not a solution of Eq. (1), exists. It may not, however, be interpreted 
as a canonical momentum. Let H{p, q) be a function in phase space which is the 
hamiltonian of the system, whenever p and q are canonical pairs. It is related to 
through the Legendre transformation. 

For a given C, Eq. (2) is an algebraic equation for H. One may, however, take 
a different point of view. For a given functional form of H{p, q), Eq. (2) may 
be considered as a differential equation for C. Its non unique solutions differ 
from one another by total time derivatives. One may also study the same system 
in the momentum space. Let jC^ {p, p) be a lagrangian in p space. It is related to 
H{p, q) as follows, 

H{p,^) = -p^+C^{p,p). (3) 

Here the functional dependence of H on its argument is the same as in Eq. (2). In 
principle, Eq. (3) should be solvable for up to an additive total time derivative 
term. Once is known the actual trajectories in p space are obtainable from an 
Euler-Lagrange equation analogous to Eq. (1) in which q is replaced by p. That 
is 

d dCP _djCP _^ 
dt dp dp 

The derivative along an actual p trajectory is the canonical coordinate 
conjugate to p . Calculated on a virtual orbit, it is not. 

A formulation of quantum statistical mechanics is possible by extending the 
conventional phase space and by applying the canonical quantization procedure 
to the extended quantities in this space. Assuming the phase space coordinates p 
and q to be independent variables on the virtual trajectories, allows one to define 
momenta np and iTq, conjugate to p and q, respectively. One may combine the 
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two pictures and define an extended lagrangian in the phase space as the sum of 
p and q lagrangians, 

^p, q,p, q) = -m -qp + ^''{p,p) + ^'^{q, <?)• (5) 

The first two terms in Eq. (5) constitute a total time derivative. The equations of 
motion are 

ddC dC d dCP dCP „ 



dt dp dp dt dp dp 
ddC dC d dC dC 



= 0, (6) 

dt dq dq dt dq dq ^' 
The p and q in Eqs. (5-7) are not, in general, canonical pairs. They are so only 
on actual trajectories and through a proper choice of the initial values. This gives 
the freedom of introducing a second set of canonical momenta for both p and q. 
One does this through the extended lagrangian. Thus 

dC dCP 

TTp = ^ = -jp- - q, (8) 

dp dp 



T^<i = — = — -P- (9) 



dC _ dCi 
dq dq 

Evidently, TTp and TTg vanish on actual trajectories and remain non zero on virtual 
ones. From these extended momenta, one defines an extended hamiltonian, 



'H{-Kp,'Kq,p,q) = pKp + qiTq - C = H{p + Trq,q) - H{p,q + 'Kp) 



y- 

^ Til 



a^-. - 5^-. ) ■ (10) 



Using the canonical quantization rule, the following postulates are outlined: a) 
Let p, g, TTp and tt, be operators in a Hilbert space, X, of all square integrable 
complex functions, satisfying the following commutation relations 

d 

K,9] = -i?i, Trg = -ih—, (11) 

d 

[TTp,p] = -ih, TTp = -iflg^, (12) 

[p,q] = [7Tp,7Tq]^0. (13) 

By virtue of Eqs. (1 1-13), the extended hamiltonian, H, will also be an operator 
in X. b) A state function x{p, q,t) € X is assumed to satisfy the following 
dynamical equation 

dy d d 

ih— = nx=[H{p-ih—,q)-H{p,q-ih—)]x 

~ ^~n\ \~dp^dq^~~dq^dp^j^' 
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c) The averaging rale for an observable 0{p, q), a c-number operator in this 
formalism, is given as 

< 0{p, q) >= j Oip, q)x*{p, q, t)dpdq. (15) 
To find the solutions for Eq. (14) one may assume 

X{p,q,t) = F{p,q,t)e-^P"/'^ (16) 

The phase factor comes out due to the total derivative in the lagrangian of Eq. 
(5), -d(pq)/dt. The effect is the appearance of a phase factor, exp{—ipq/h), in 
the state function that would have been in the absence of the total derivative. It 
is easily verified that 

(p-i?i^)x = *^^e-^f^/'^ (17) 

[q-in—)x = ih—e-'^'"' (18) 

Substituting Eqs. (17) and (18) in Eq. (14) and eliminating the exponential 
factor gives 

3 3 3F 

H{-ih—,q)-Hip,-ih—)F = ih—. (19) 
3q 3p at 

Equation (19) has separable solutions of the from 

F(p,q,t)=ij{q,t)r{p,t), (20) 

where ip{q, t) and <ji)(p, t) are the solutions of the Schrodinger equation in q and 
p representations, respectively. The solution of the form (16) associated with 
anti-standard ordering rale satisfies Eq. (14) and is one possible distribution 
function. For more details on the admissibility of the distribution functions, their 
interesting properties and the correspondence rales, one may consult Sobouti 
and Nasiri [2]. 

3 Generalized distribution function and ordering rule 

In 1932 Wigner [1] proposed the distribution 

= i / y^p^^i + y> ^'"'^''dy' (21) 

for a system in mixed state represented by a density matrix p. The expectation 
value for an operator O calculated with W{q, p) has the same value as ordinary 
quantum average with wave function -0, i.e. 

<6>^= Jw{q,p)Oiq,p)dpdq, (22) 
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where 0(q,p) is a classical function corresponding to operator O and is given 
according to wigner prescription by[3] 

Oiq,p) =J<1- + I > e'P^/'^dz. (23) 

Similarly, the distribution function of Eq. (16) could be written as 

xiq,P) = ^ / < 1\P\1 + ^ > e^P^'^dz, (24) 
and corresponding averaging rule 

< 6{q,p) >x= j x{q,P)0{q,p)dpdq, (25) 
where the classical function 0{q,p) is given by 

0{q,p) = j < q\0\q + z> e'P'/^dz. (26) 
Let a generaUzed distribution function be defined as 

Pa{q,p) = ^J <1 + (^Af>\q + (a + 1)^ > e'^^'^dz, (27) 
and the corresponding ordering rule as 

Oa{q,p) = J <q + az\d\q+{a + l)z > e'^'/'^dz, (28) 

where a is a parameter specifying the given distribution function, that is, for 
a = ^ one has the Wigner function, for a = one has the standard distribution 
function and etc. The origin of a is in the canonical transformations in the 
extended phase space that serve to obtain the different distribution functions [2]. 
Equations (27)and (28) may be rewritten as 



Pa{q,p) = ^ / < q\e'^"^/''pe-'P°'^/''\q + z> e'^'/^dz 
Oa{q,p) = j < qW^°''''^de-'P°'''^\q + z > e'P'/^dz. 



(29) 



(30) 



Equations (29) and (30) relate the generalized distribution function and the cor- 
responding ordering rule to those of the (SN) representation by a unitary trans- 
formation. Thus, using the EPS representation as the basis, one may introduce 
different representations using Eqs. (29) and (30) for different values of a. As 
an example, consider an operator as O = q"^p"'. Using Eq. (28) the classical 
function corresponding to this operator after a lengthy calculations becomes 

Oa{p,q) = ^^^T(iM>""'^a™-^n < r. (31) 

\n — r)\ 

Assuming a = ^ , Eq. (31) will be identical with that of Wyle- Wigner ordering 
rule [3], and for a = 0, one obtains the standard ordering rule [4]. 
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4 Conclusions 

In the phase space approach to quantum mechanics each distribution function 
corresponds to an ordering rule. In this paper we introduce a generalized dis- 
tribution function and corresponding generalized ordering rule. A parameter, a 
serves to specify the different quantum distribution functions. The origin of the 
parameter a comes from the extended canonical transformations in the extended 
phase space. The cases of SN and Wigner distribution functions are worked out 
as examples. 
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